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ABSTRACT 

Leo P is a low-luminosity dwarf galaxy discovered through the blind Hi Arecibo Legacy Fast ALFA 
(ALFALFA) survey. The Hi and follow-up optical observations have shown that Leo P is a gas-rich 
dwarf galaxy with active star formation, an underlying older population, and an extremely low oxygen 
abundance. We have obtained optical imaging with the Hubble Space Telescope to two magnitudes 
below the red clump in order to study the evolution of Leo P. We refine the distance measurement 
to Leo P to be 1.62 ± 0.15 Mpc, based on the luminosity of the horizontal branch stars and 10 newly 
identified RR Lyrae candidates. This places the galaxy at the edge of the Local Group, ^ 0.4 Mpc 
from Sextans B, the nearest galaxy in the NGC 3109 association of dwarf galaxies of which Leo P 
is clearly a member. The star responsible for ionizing the Hii region is most likely an 07V or 08V 
spectral type, with a stellar mass > 25 Mq . The presence of this star provides observational evidence 
that massive stars at the upper-end of the initial mass function are capable of being formed at star 
formation rates as low as ^ 10“® Mq yr“^. The best-fitting star formation history derived from the 
resolved stellar populations of Leo P using the latest PARSEC models shows a relatively constant star 
formation rate over the lifetime of the galaxy. The modeled luminosity characteristics of Leo P at early 
times are consistent with low-luminosity dSph Milky Way satellites, suggesting that Leo P is what a 
low-mass dSph would look like if it evolved in isolation and retained its gas. Despite the very low mass 
of Leo P, the imprint of reionization on its star formation history is subtle at best, and consistent with 
being totally negligible. The isolation of Leo P, and the total quenching of star formation of Milky 
Way satellites of similar mass, implies that local environment dominates the quenching of the Milky 
Way satellites. 

Subject headings: galaxies: dwarf - galaxies: distances and redshifts - galaxies: photometry - galaxies: 
stellar content - galaxies: fundamental parameters - galaxies: evolution 


1. INTRODUCTION 

I.l. Evolution of Very Low-Mass Galaxies 

Galaxies at the faint-end of the luminosity function 
(LF) have provided vital tests for our understanding of 
structure formation and evolution. The dearth of both 
observable lo w-mass dwarf galaxies (“missing satellite 
problem”; e.g., Kauffmann et al.|1993[|Klypin ^ al.|1999 ' 
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Moore et al.|1999 1 and of hi gher surface brightness dwarf 
galaxies (“too big to fail”; Boylan-Kolchin et al. 2011| 
compared to simulations based on the ACDM model 
has proven to be a multi-decade challenge. These issues 
have been partly alleviated with higher resolutio n simu- 


et al. 


lations that include more baryonic effects (e.g., Brooks 
IBenitez-Llambay et al.|[2014| |Sawala et al. 


2015 


Ohorbe et al.|2015). In addition, large sky surveys 


have detected an increasing number of satellite systems 
providing progressively better agreement between theory 
and observations for the missing satellite problem (e.g., 


Koposov et al. 20151. While these results have reduced 
the inconsistencies, future revisions to structure forma¬ 
tion theories will undoubtedly be anchored by observa¬ 
tional constraints from studies of very low-mass systems. 

Furthermore, galaxies at the faint-end of the LF are 
broadening our understanding of star formation and 
galaxy evolution. The star formation properties of gas- 
rich, low-mass galaxies encompass a larger range of pa- 


“main-sequence” galaxies (Brinchmann et al.|2004 Both- 

well et al.||2009[ 

Huang et al.| 

p012| IMcQuinn et al. 

2015b|). It is unc 

ear whether t 

tiis is due to an inher- 


gas column densities and metallicities, or to a variation 
in the primary processes that regulate star formation. 

In the nearby Universe, much work has been done to 
understand the evolutionary history of the rich popula- 


tion of nearby low-mass galaxies (e.g., Skillman et al. 
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Figure 1. Three-color HST image of Leo P showing both the main star forming complex and underlying older stellar population extending 
to greater radii. The image was created by combining the F475W band image (Blue), the average of F475W and F814W band images 
(Green), and the F814W band image (Red). The field of view shown here is approximately 3.1'xl.8'. Note the dither pattern used did 
not cover the chip gap which is visible across the top of the image. 

metallicity (3% solar), significant gas content (8.1 x 10^ 
Mq), proximity (1.62 Mpc), and isolation (0.4 Mpc from 
its nearest neighbor and outside of any group environ¬ 
ment) present a unique test for our understanding of both 
structure formation and galaxy evolution. Comparable 
in luminosity to some of the dSph satellites of the Milky 
Way and a few of the dirrs inside the LG, Leo P pro¬ 
vides an opportune target for studying evolution at the 
faint end of the LF in a simplified system. Not only are 
many complex variables eliminated (e.g., tidal and ram 
pressure stripping, significant possible merger events, etc. 
that are present in the majority of galaxies detected in 
this luminosity regime), but Leo P crosses over into the 
mass regime where (i) theories of structure formation 
predict difficulty in galaxies retaini ng their baryons due 


cause relative differences between low-mass galaxies can 
depend significantly on large-scale environmental dif¬ 
ferences. This latter factor is often described as the 
morphology-density relation, or the prevalence of gas- 
poor dwarfs (i.e., dwarf spheroidals or dSphs) to be lo¬ 
cated in higher galactic density environments whereas 
gas-rich dwarfs (i.e., dwarf irregulars or dirrs) are pre- 
domi nantly found in lower density field environments 


(e-K-, 

Einasto et al. 

[1974 

Binggeli 

1987 

van den Bergh 

1994 

Cote et al.|20(J 

9The closest (and most accessible) 


and thus disentangling the impact of the internal versus 
external evolutionary processes is a significant challenge. 
Ideally, one could study the properties of a low-mass sys¬ 
tem located outside of a group environment. While such 
low-mass, low-luminosity systems are expected to be nu¬ 
merous, they have generally eluded detection outside of 
the immediate vicinity in the Local Group (LG) due to 
their faint and sometimes low surface brightness nature. 

An exception is Leo P, a low-luminosity (My = —9.27 
mag), gas-rich galaxy discovered just outside the zero- 
velocity boundary of the LG. Leo P’s combined proper¬ 
ties of extremely low stellar mass (5.6 x 10^ M©) and 


to stellar feedback processes (e.g. , Larson 1974 Dekel[ 
fc Silk||1986[ |Mac Low &: Ferrara 1999| Ferrara fe Tol-| 
stoy|2000p and (ii) theories of the impact of reionization 


predict the photo-evaporization of the gas content and a 
sudden and irrevocable quenching of star formation (e.g.. 
Babul & Rees|1992[|Efstathiou||1992||Thoul & Weinberg] 

1996||. 

In this work, we present a detailed study of the star 
formation history (SFH) and evolution of Leo P using 
HST optical imaging of Leo P that reaches ~2 mag be¬ 
low the red clump (RC). In Section 2 we discuss the 
observations, data reduction, and variable star identifi¬ 
cation. In Section 3, we refine the distance measurement 
to Leo P using analysis of horizontal branch (HB) stellar 
populations and variable stars. In Section 4, we describe 
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the methodology for deriving the best-fitting SFH, and 
in Section 5 we present measurements of the SFH and 
chemical evolution of Leo P. In Section 6 we present a 
comparison of the evolution of Leo P with previously 
studied Milky Way dSph galaxy satellites and three LG 
dirrs. Our conclusions are summarized in Section 7. 

1.2. A Brief Synopsis of Leo P 


Leo P (AGC 208583; Giovanelli et al. 20131 was dis- 

EW survey 


(ALFALFA; 

Giovanelli et al. 

2005 

Haynes et al. 

2011) 

The ALFALl 

hA survey is a b 

ind extragalactic survey ii 


the Hi 21 cm line covering over 7000 square degrees of 
high Galactic latitude sky. Follow-up optical observa¬ 
tions at the location of the ALFALFA Hi detection, in¬ 
cluding WIYN 3.5m BVR and KPNO 2.1m Ha imaging, 
confirmed the presence of bo th a resolved stellar popula¬ 
tion and a single Hii region ( Rhode et al.||2013 |. 

Optical spectroscopy of the Hii region enabled a direct 
measurement of the auroral [O HI] A4363 line, yield ing an 
oxygen abundan ce of 12 -|- log(0/H) = 7.17±0.04 (| Skill- 
man et ar]|2013), ^ 3% Zq (based on a solar abu ndance 


of 12 -I- log(U/H) = 8.68; Asplund et al.||2009|). This 
abundance measurement showed that Leo r is the low¬ 
est metallicity, gas-rich galaxy known in the Local Vol¬ 
ume. Its properties are c onsistent with the luminosity- 
metallicity relationship in Berg et al.| (2012|). 

Deeper optical imaging of the resolved stellar popula¬ 
tions from the Large Binocular Telescope (LBT) enabled 
a distance measurement based on the tip of the red giant 
branch (TRGB) detection of 1.72lg;4g (McQuinn et al. 
20131. The large lower uncertainty is due to the sparse- 
ness of the red giant branch (RGB) stars in such a low- 
mass galaxy. This distance measurement placed Leo P 
just outside the zero velocity boundary of the LG. 

Hi spectral line imaging from the VLA shows a small 
amplitude rotation (14 = 15 ± 5 km s“^), with no obvi¬ 
ous sig ns of interaction or in-falling ga s at larger spatial 
scales (Bernstein-Gooper et al. 2014|. These observa¬ 
tions revealed that Leo P has one of the lowest neutral 
gas mass of any known low-metallicity dwarf, with a mass 
ratio of gas to stars of 2 : 1 and total mass to baryonic 
mass of > 15 : 1. Follow-up high-sensitivity observations 
with GARMA of the GO (J= 1 —0) transition did not 
detect any GO emission, but placed stringent upper lim¬ 
its on the CO luminosity of L^o < 2930 K km s“^ pc^ 
(Warren et al. submitted). The molecular hydrogen mass 
remains uncertain as the ratio of CO to H 2 is not well 
constrained at low-metallicities. In Table [T] we summa¬ 
rize the basic properties measured in Leo P from these 
previous studies. 

2. DATA PROCESSING 
2.1. Observations 

HST observations of Leo P were obtained usin g the 
Advance C amera for Surveys (ACS) instrument (Ford 
et al.|1998 ) in the F475W (Sloan g) and F814W (I) filters 
between April 23 and 26 2014 (GO 13376; PI McQuinn). 
The observations were taken over 17 orbits grouped in 
visits of 2 — 3 orbits, and followed the ultra-deep field 
dither pattern between exposures, shifted by 2 — 3 pixels 
to aid in the removal of hot pixels and to average-out 
the detector response. The observations were designed 


to reach a photometric depth of ~ 2 mag below the RC 
to aid in (i) constraining the SFH at older times and (ii) 
measuring the distance from the HB stars and RR Lyrae 
stars. 

Similar to the successf ul Local Cosmology from Iso- 


Skillman et al. 


20141 program. 


lated Dwarf^(LCID, e.g. 
our observation strategy was designed to optimally sam- 
ple the light curves of sh ort period variable stars (e.g., 
Bernard et al.|2008 20091, such as RR Lyrae, which pro- 
vide independent constraints on the distance to Leo P 
and the SFH. To achieve this, each orbit was split be¬ 
tween the two filters, with an alternating pattern for the 
exposures (i.e., F475W, F814W, F814W, F475W, etc.). 
In addition, the visits were executed sequentially and 
grouped within a 3-day time period, while minimizing 
orbits with a ^ 12 hour cadence. These observations 
thus maximize sampling within the 17 orbits of a ~ 0.6 
day period RR Lyrae star across 6 periods with minimal 
repetition of sampling at the same point in a light curve. 
The observational details are summarized in Table [TJ 

HST Drizzlepac v2.0 was used to create mosaics in 
each filter. First, the charge transfer efficiency corrected 
images (i.e., flc files) were processed with Astrodriz- 
ZLE to remove cosmic-rays (CR), creating a set of CR 
cleaned images. Second, these CR cleaned images were 
processed through the task tweakreg to measure the 
astrometric shifts between the images. Third, the orig¬ 
inal flc.fits images were processed with Astrodriz- 
ZLE using the tweakreg shifts to create a combined 
mosaic in each filter. 

In Figure we present a color image of Leo P from 
the HST AGs observations created using the F475W and 
F814W mosaics, along with an averaged image of the two 
filters. As seen in Figure [T] the galaxy has an irregular 
stellar morphology that is well-resolved and a single Hii 
region clearly identifiable in the image. 

2.2. Photometry and Color-Magnitude Diagram 

Point spread function (PSF) photometry was per¬ 
formed on the flat, charge transfer efficiency corrected 
images (flc.fits) with the DDLPHOT photometry pack¬ 
age, a m odified version of HSTPHOT with an ACS specific 
module (Dolphin 2000). The photometric catalog was 
filtered for well-recovered stars (i.e., output error flag 
< 8) and with a signal-to-noise ratio > 5. Point sources 
with high sharpness or crowding values were rejected 
(i.e., iy sharp~tAsharpf' > 0.075; (V crow crowd )> 0.8). 
Sharpness indicates whether a point source is too broad 
(such as background galaxies) or too sharp (such as cos¬ 
mic rays). Crowding measures how much brighter a star 
would be if nearby stars had not been fitted simultane¬ 
ously; stars with higher values of crowding have higher 
photometric uncertainties. Spatial cuts were applied to 
the filtered photometric catalog. These cuts were based 
on the ellipticity (e =0.52), semi-major axis (1'.2), and 
posit ion angle (PA = 335°) of Leo P (McQuinn et al. 
|2013| . The final elliptical parameters are listed in 'Fa¬ 
ble [l] We expanded the spatial cuts for variable star 
identification (see Section 2.4), extending the search to 
include the halo of the galaxy. 

Artificial star tests were performed to measure the 
completeness limits of the images using the same pho- 

^ http://www.iac.es/project/LCID/ 
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Table 1 

Leo P Properties and Observations 


Parameter 

Value 

Updated Distance- 
Dependent Values 

R.A.(J2000) 

10:21:45.1 


Decl. (J2000) 

+18:05:17 


Distance (Mpc) 

79-^0.14 

-^•'^-0.40 

1.62 + 0.15 

Mv (mag) 

—Q 41+0.17 

-9.27 + 0.20 

M* (Mq) from M/L 

5.7l?-^ X lOS 


Af, (M©) from SFH 


5.61°'^ X 10® 

M*/Ly 


1.25 

(Mq) 

9.3 X 10® 

8.1 X 10® 

Vsys (km s“l) 

260.8 + 2.5 


Leo (K km pc^) 

< 3300 

< 2930 

12-flog(0/H) 

7.17+0.04 


Lhc (erg s-Q 

6.2 X 10®® 

5.5 X 10®® 

SFRhc (Mq yr-l) 

4.9 X 10-® 

4.3 X 10“® 

Av (mag) 

0.071 


Semi-major axis (') 

1.2 


Ellipticity 

0.52 


Position angle (°) 

335 


HST Program ID 

HST-GO-13376 


AGS Filters 

F475W; F814W 


Exposure F475W (s) 

25,829 


Exposure F814W (s) 

18,496 



Note. — Sum mary of the properties o f Leo P based on measure¬ 


ments reported in 


G tovaneUietab] 2013t; [Rhode ^ aI.ll|201 3| jM^ 


I Quinn et al.lll201al);lakillman et al. f! 
||:^U14[| and Warreri et al. isubmittedj 


2Uiii|l ; ItSernstem-Cooper et al. 
and the revised distance mea- 


is in the Local Standard of Rest Kinematic frame. Ho-based SFR 

is based on the calibration by|Kennicutt et al.|( 

|1994|. Foreground 

extinction estimate is based on the dust maos i 

Tom ISchlegel et al.l 


||1998[l with recalibration from|Schlafly fc Finkbein'S^ (|2UllTr 


tometry package. Approximately 600k artificial stars 
were injected into the individual images following the 
spatial distribution of all sources in the non-Hltered pho¬ 
tometric catalog. The final artificial star lists were fil¬ 
tered using the same parameters that were used to pro¬ 
duce the final photometry catalog. 

Figure shows the CMD from the final photometry 
plotted to the 50% completeness level as determined from 
the artihcial star tests. Representative uncertainties per 
magnitude are plotted and include uncertainties from the 
photometry and uncertainties recovered from the artih¬ 
cial star tests. The photometry was corrected for ex¬ 
tinction (Ai? 475 w = 0.086 mag; Afsmw = 0.039 mag) 
based on the dust maps of ISchlegel et ah ] (|1998 |) with 
recalibration from Schlahy fe Finkbeiner| (|2011|). The 


main sequence (MS), RGR, HR, and the RC are identi- 
hable in the CMD. In addition, there is a sparse popula¬ 
tion of red and blue helium burning stars (i.e., blue loop 
stars with F814W < 25 mag). These stars are intermedi¬ 
ate mass stars (2 M 0 < M <15 Mq) and are unambigu- 
ous signs of recent (t< 500 Myr) star-formation activity 
(Dohm-Palmer & Skillman 2002 McQuinn et al. 2011). 
We investigate the star-forming properties in Leo F in 
Section 5. 

Figure shows the spatial distribution of stars sepa¬ 
rated into course age bins of “young” and “old” stars. 
Stars were classihed as young and old based on their po¬ 
sition in the CMD. Specihcally, young stars were selected 
from the upper MS corresponding to a B or O spectral- 
type classihcation with lifetimes < 300 Myr (i.e., F475W 



Figure 2. The CMD of Leo P from HST ACS imaging plotted to 
the 50% completeness limit determined from artificial st ar tests. 
The photometry w as corrected for Galactic extinction (|Schlafly| 
I&: Finkbeiner||2011|l. Representative uncertainties per magnitude 
bin are plotted and include both photometric and completeness 
uncertainties. The MS, RGB, HB, and RC are seen in the CMD. 

< 25.7 mag and F475W — F814W < 0.2 mag), and the 
blue and red helium burning branch stars which have 
a similar maximum lifetime. In Figure we overlay 
an isochrone of 300 Myr to demonstrate the location of 
young stars in the CMD. The old category is populated 
with RGB stars identified from the CMD. Despite the 
coarse time bins, the selection of young vs. old stars is ro¬ 
bust. In Figurej^ the young stars are more centrally con¬ 
centrated with an irregular distribution. A small number 
of young star candidates are located in the outer regions 
of the galaxy. The locations of the RR Lyrae candidates, 
identified from the photometry, are also shown (see Sec¬ 
tion 2.4). 


2.3. Resolved Stars in the Hn Region 

The area within the one Hii region in the galaxy is re¬ 
solved into 7 individual point sources in the HST images. 
These sources were blended in previous, ground-based 
imaging of Leo P, leading to some speculation on the na¬ 
ture of t he star(s) responsib le for ionizing the surround¬ 


ing gas (Rhode et al. 2013 see their §4.1.2). Figure]^ 


shows the CiVID of Leo P with the F475W magnitude as 
the ordinate axis. Here, we highlight the stars located at 
the position of the Hii region in Leo P in red, including 
the luminous star likely responsible for the ionization of 
the region. Note that these stars are bluer than expected 
in the CMD, which is most likely due to contamination 
by emission lines in the Hii region. The most luminous 
of these stars has de-reddened F475W and F8I4W mag 
of 21.509 ± 0.001 and 21.892 ± 0.002 respectively. Us¬ 
ing the calibration of Saha et al. (2011) and an adopted 
distance modulus of 26.05 ± 0.20 (see Section 2.4 for dis¬ 
cussion of the calibration and Section 3.3 for the adopted 
distance), the ACS filter luminosities transform to an ab- 
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Figure 3. Spatial distribution of “young” (age < 300 Myr) and 
“old” (age > 300 Myr) stars in Leo P as defined in the text. The 
young stars are more centrally concentrated, with an irregular dis¬ 
tribution. Despite the central concentration, there are still young 
stars located in the outer regions of Leo P. 


solute Johnson V mag of —4.43. A star of this magnitude 
in the upper main sequence region of the CMD is con¬ 
sistent with an 05V spectral type, according to the em- 
pirical catalog of absolute magnitudes of OB stars from 
Wegner I (|2000l see their Table 1 of smoothed absolute 
magnitudes). However, if this star is truly an 05V star, 
then the star with a similar color but ^ 1 mag fainter 
at F475W = 23.03 mag would correspond to a star with 
an 09.5 spectral type, which would be hot enough to 
create a second Hii region. Since no Hii region is found 
around this fainter star, it is likely that both of these 
point sources are blended binary stars. High mass stars 
have been shown to have b inary fractions in excess of 
50% ( Kobulnicky et al.|2014 ). If we assume each of these 
MS stars is comprised of equal mass binaries, the bright¬ 
est source located in the known Hii region could have 
two constituent 07 or 08 stars. The fainter source con¬ 
stituents would then be B2 or B3 stars, just below the 
li mit for hosting a de tectable Hii region. 

Rhode et al.j (2013) calculated a SFR based on the Ha 
luminosity of the Hii region of 4.9 x 10“® Mq yr“^. Us¬ 
ing our adopted distance to Leo P, this SFR is revised 
downward to a value of 4.3 x 10“^ Mq yr“^, compared 
with our SFR derived from the stellar populations over 
the past 4 Myr of 2.1 x 10“® Mq yr“^ (see Section 5). 
In this regime, the Ha luminosity has been shown to be 
a poor tracer of the actual star-formation activity due, 
in part, to the incomplete sampling of the high mass end 


of the initial mass function (IMF) (e.g., Boselli et al. 

2009 Goddard et al. 

20101 Koda et al. 20121 

Weidner 

& Kroupa I 

2005) and Phamm-Altenburg et a 

nWT 


have suggested there may be an upper limit to the max¬ 
imum stellar mass that can form at low SFRs. In this 
scenario, the slope of the upper end of the I MF varies as 
a function of S FR. Using the framework of|Weidner fc 


Kroupa 


et al. 


(2005 their Figure 11), and Phamm-Altenburg 


(2()07[ their Figure 4), the maximum stellar mass 


that would be able to form at a SFR = 10“^ Mq yr“^ is 
^ 2.5 Mq. This is an order of magnitude lower than the 
mass of an ~08 star (M* = 25 Mq) likely responsible 
for ionizing the Hii region. Indeed, the mere presence of 


the Hii region in Leo P provides observational evidence 
that the maximum stellar mass capable of being formed 
at low SFRs is significantly greater than 2.5 Mq. 



F475W-F814W (mag) 


Figure 4. CMD with the F475W plotted as the ordinate axis. 
Stars located at the position of the Hll region are circled in red. 
PARSEC stellar evolution isochrones are over plotted for a 12 Gyr 
population with metallicities varying between [Fe/H] of —1.9 to 
— 1.7 and for a 300 Myr population with a metallicity of [Fe/H] 
= —1.5, as shown in the legend. 


2.4. Variable Star Identification 

We used the multi-epoch photometry to search for 
short-period variable stars. As noted in Section 2.2, the 
spatial cuts were relaxed to include a larger search area 
as RR Lyrae variable stars can be located in the stellar 
haloes of galaxies. Variable star candidates were identi¬ 
fied first using a series of auto mated cuts, similar t o those 
used in the study of IC 1613 ( Dolphin et al.|2(iM |. Each 
variable star was required to meet hve criteria testing 
the photometry, variability, and periodicity. First, the 
object had to have well-measured photometery. While 
it is possible for a variable star to be part of a blend, 
our PSF-fitting photometry attempts to fit the profile to 
that of a single star and thus is unreliable for measuring 
blended stars. The next cuts were intended to eliminate 
non-variable or weakly-variable stars. The simplest re¬ 
quirement was that the RMS scatter of the magnitude 
measurements, 


^mag — ( 


N 


F475W 


1 

+ -^F814VF 


Nf475W 

*( {F475W^-F475Wf 

2=1 


NfSIAW 

+ {FSUWi - FSUWf))^^^ 

2=1 


(1) 
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had to be 0.14 magnitudes or greater. Next, the overall 
reduced 


X^ = 


1 


Nf475W + Nfs14W 


Nf475W 

( E 


{F475Wi - FA75Wy 


N'f814W 


^ {F8UW, - F8UW)^ ^ 


had to be 3.0 or greater. These two cuts eliminate 
stars that are weakly variable and those for which the 
photometry signal-to-noise is insufficient to discern vari¬ 
ability. We also clipped 1/3 of the extreme points and 
required that the recalculated reduced exceed 0.5; this 
was intended to eliminate non-variable stars with a few 
bad points. 

The final test was a m odified Lafler-Kinman algorithm 
(Lafler & Kinman 1965), which tests for periodicity. This 
was implemented by computing 0 for periods between 
0.1 and 4.0 days. The 0 parameter is calculated by de¬ 
termining the light curve for a trial period and using the 
equation 




0 = 




^^2 


( 3 ) 


where N is the number of exposures for a given filter, 
rrii is the magnitude at phase i, and m is the mean mag¬ 
nitude. If the trial period is the correct period, each 
magnitude rrii will be close to the adjacent magnitude 
rrii+i, giving a value of 0 that decreases as 1/N^. If the 
trial period is incorrect, there will be less correlation and 
consequently larger values of 0. Because we had data in 
two filters, we combined the 0 values with: 



0 .6 1 1.5 2 0 .6 1 1.5 2 

Phase Phase 


Figure 5. Light curves for the ten identified RR Lyrae candidate 
stars in Leo P. The measured periods are noted in the individual 
panels for each star. Complete photometry is available in Table 


Table 2 

Photometry of RR Lyrae Candidates 


0 = 


(i/ye 


F475W 


i/y© 


F814W 


( 4 ) 


For a star to pass the periodicity test, its minimum value 
of 0 had to be 0.85 or less. For more discussion on the 
use of the modified Lafler-Kinman algorithm t o test for 
periodicity and sp ecifics on the parameters, see |Dolphin 
et al.| (1200211200^. 

Uur photometry identified 13448 total objects, of 
which 45 passed our automatic selection process. Three 
of these were clearly RR Lyrae stars by location in the 
CMD, light curve shape, and period. The criterion 
was relaxed, and the search was constrained to an area 
centered on the horizontal branch, and ten more candi¬ 
dates were identified. Of these new ten, three of these 
candidate variables appeared to not be clean detections 
upon manual examination of the images and were re¬ 
moved from the list resulting in a final total of ten candi¬ 
date RR Lyrae stars. The light curves for the ten variable 
stars are presented in Figure complete photometry in 
the ACS filters is available in Table [21 
The absolute magnitudes of RR Lyrae stars have been 
determined based on the Johnson BVRI hlter systems. 
Therefore, we transform our measured ACS F475W and 
F8I4W magnitudes to standard J ohnson V and I m ag- 
nitudes using the calibration from Saha et al. (2011): 


Star ID 1 

10:21:44.388 +18:05:07.63 

Epoch 

F475W 

(mag) 

56770.631305 

56770.710315 

56770.829933 

56770.909546 

56771.626484 

26.70 ± 0.07 
27.08 ± 0.08 
27.12 ± 0.10 
27.26 ± 0.10 
27.04 ± 0.09 

Epoch 

F814W 

(mag) 

56770.647902 

56770.692595 

56770.846531 

56770.891826 

56771.643081 

26.04 ± 0.10 
25.95 ± 0.10 
26.03 ± 0.09 
26.48 ± 0.15 
26.25 ± 0.11 


Note. — Sample photometry of the first five epochs for Star 
ID 1 in both ACS filters. The RA and Dec (J2000) coordinates 
for the star are listed at the top. These magnitudes have not 
been corrected for foreground extinction. The full seventeen epoch 
photometry for all ten RR Lyrae candidates is published in its 
entirety in the electronic edition of this journal. A portion is shown 
here for guidance regarding its content. Mean magnitudes for each 
of the ten stars can be found in Table 1^ 
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V = 0.026 -h F4751U - 0.406 * (F475W - F8UW) (5) 
I = -0.038+ F814W +0.014 *(F475W-F814W). (6) 


These transformations are an improvement over the pre¬ 
vious calibrations from Sirianni et al. (20051 as they were 
derived based on a larger number ot stars down to fainter 
magnitudes and over a greater range in color. Fu rther- 
more, the transformations from Saha et al. (2011) were 
based on updated corrections for time-dependent charge 
transfer losses in the ACS detector. 

To compute mean magnitudes, we calculated a phase- 
weighted average using 


N 


r) =-2.51og^ii±—(7) 


where (j)i is the phase and rrii is the magnitude at each 
point along the light curve. Finally, the mean magni¬ 
tudes for all ten stars are computed to be F475W = 
26.88±0.02, F814W = 26.18±0.07, and V = 26.62±0.03. 
Note, these mean magnitudes are not corrected for fore¬ 
ground extinction. Table summarizes the mean mag¬ 
nitudes for each star in the three filters and the final 
mean magnitudes for all ten stars. The location of the 
RR Lyrae candidates within Leo P are shown in Figurej^ 



F475W-F814W (mag) 


Figure 6. A magnified view of the HB and RC region of the 
CMD. A portion of the RGB is identifiable to the right. The ten 
RR Lyrae candidates are plotted in red. As expected, these stars 
lie in a tight sequence in CMD space. 


Figure^shows the location of the ten RR Lyrae candi¬ 
date variable stars in the CMD. All ten stars were found 
in a narrow range in luminosity. Because the HB is ^ 
2 magnitudes above the 50% completeness limit for our 
observations, we likely did not miss any RR Lyrae stars 
due to depth of photometry. Also, because of the number 
and spacing of the epochs, we are likely to have a very 
high completeness. The data provided excellent coverage 


for the range of relevant variable star periods. Due to the 
large number of observational epochs and the fact that 
the cadence was set to minimize redundancies, the max¬ 
imum gap in the phase coverage was always less than 
0.15 dws for periods in the range 0.1 to 1.0 days (see 
Figure!^. Consequently, our mean magnitude measure¬ 
ments are robust, as indicated by the very small disper¬ 
sion [<7v = ±0.03) in the sample. Thus, the distance to 
Leo P derived from the RR Lyrae stars (see Section 3) is 
limited by the systematic uncertainty of the calibration 
used. 


3. THE DISTANCE TO LEO P 

The previously published distance modulus mea sure- 
ment of 26.19^Q gg mag from McQuinn et al. (2013) was 
determined by applying the IRCB distance method to 
the photometry of the resolved stellar population from 
LBT V and I band imaging. Despite photometry reach¬ 
ing 3 magnitudes below the TRGB, the distance uncer¬ 
tainties in the ground-based observations are large due 
to the small number of stars in the RGB. Using Monte 
Carlo simulations of synt hetic stellar population s analo¬ 
gous to low-mass galaxies, |McQuinn et al.|(|2013|) showed 
that the measured break in the LBT 1 band LF at the 
top of the RGB may, in fact, be below the actual TRGB 
luminosity. The large lower uncertainty estimated from 
the simulations put constraints on the actual TRGB of 
up to 0.5 mag brighter than the measured value. Us¬ 
ing the same TRG B distance method described in |Mc-| 
Quinn et al. (2013), we re-measured the TRGB distance 
to Leo P using the HST data. This exercise yielded an 
identical measurement of the distance modulus of 26.19 
mag as the LBT data, but does not mitigate the uncer¬ 
tainties as measured by the Monte Carlo simulations. 

However, the new HST data reach 2 mag below the 
HB, enabling two additional distance techniques to be 
applied. First, as described above, the cadence of our 
observations was designed to maximize identification of 
RR Lyrae stars, thus enabling a distance measurement 
based on the average V-band luminosity of this class of 
stars. Second, the luminosity of the HB provides a means 
to measure the distance based on a standard candle ap¬ 
proach. The HB luminosity can be more robustly mea¬ 
sured than the TRGB with sparse data as the mean lu¬ 
minosity of the HB is not as impacted by the presence 
or absence of a few individual stars. 


3.1. RR Lyrae Distance Measurement 

As discussed above, the mean V magnitude of the ten 
RR Lyrae stars is 26.62 ± 0.03 mag. To account for fore¬ 
ground extinction, w e subtract Ay = 0.07 1 mag based 
on the dust ma ps of Schlegel et al. (19 98) w ith the re¬ 
calibration from Schlahy & Finkbeiner iSoil I. 

There has been signihcant work done to calibrate the 
luminosity of the RR Lyrae stars as a function of metal- 
licity, however the zero-point can still vary by ~ 0.1 


mag between calibrations (for detailed reviews, see Smith 
[T^|Gacciari|2003[ [S andage|2Q06| [Catelan|20Q9[ 

others). Some of the differences are due to the di 


among 
ifhculty 

quantifying the luminosity dependence on stellar metal- 
licity. Regardless, the choice of calibration introduces 
a systematic error to our distance measurement, dom- 
inating the uncertaint ies. We use the calibration from 
Carretta et al. (2000) which was derived using dozens 
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Table 3 

Summary of RR Lyrae Candidates in Leo P 


Star ID 

RA 

(J2000) 

Dec 

(J2000) 

<F475W> 

(mag) 

<F814W> 

(mag) 

<V> 

(mag) 

Period 

(days) 

01 

10:21:44.388 

-1-18:05:07.63 

26.87 ± 0.09 

26.04 ± 0.07 

26.56 ±0.14 

0.61±0.06 

02 

10:21:47.080 

-1-18:03:59.16 

26.85 ± 0.08 

26.22 ±0.07 

26.62 ±0.12 

0.52±0.06 

03 

10:21:44.324 

-1-18:05:33.73 

26.89 ± 0.06 

26.11 ± 0.07 

26.56 ±0.09 

0.69±0.11 

04 

10:21:43.240 

-1-18:05:53.13 

26.90 ± 0.08 

26.16 ±0.06 

26.62 ±0.12 

0.64±0.08 

05 

10:21:44.722 

-1-18:05:39.91 

26.90 ± 0.10 

26.25 ±0.07 

26.66 ±0.15 

0.57±0.13 

06 

10:21:42.750 

-1-18:05:46.47 

26.89 ± 0.09 

26.17 ± 0.05 

26.63 ±0.13 

0.59±0.06 

07 

10:21:43.704 

-1-18:05:34.37 

26.89 ± 0.10 

26.25 ±0.07 

26.66 ±0.15 

0.61±0.04 

08 

10:21:45.227 

4-18:05:46.86 

26.91 ± 0.07 

26.11 ± 0.06 

26.61 ± 0.11 

0.57±0.04 

09 

10:21:44.633 

4-18:05:21.31 

26.84 ± 0.12 

26.27 ± 0.06 

26.64 ±0.18 

0.53±0.07 

10 

10:21:44.087 

4-18:06:22.92 

26.86 ± 0.06 

26.24 ± 0.07 

26.63 ±0.10 

0.34±0.01 


Mean Magnitudes 

26.88 ± 0.02 

26.18 ±0.07 

26.62 ±0.03 



Note. — 10 RR Lyrae variable star candidates were identified as described in Section 2.4. The coordinates, phase-weighted mean 
magnitudes, and periods of the individual stars are listed, followed by a mean magnitude computed for all ten stars. These magnitudes 
have not been corrected for foreground extinction. 


of RR Lyrae stars in globular clusters and anchored to 
several independent distance measurements: 


Mv = (0.18 ± 0.09) *([Fe/H] +1.5)+ (0.57±0.07). (8) 

As seen in Equation!^ the RR Lyrae luminosity depends 
on the metallicity oithe population. For an estimate of 


stellar evolution isochrones (PARSEC: 

Bressan et al. 

2012 

Chen et al. 

2014 

Tang et al. 

2014 

) to the RGB. 


sion of the previous Padova stellar evolution code with 
the most up to date input physics. The PARSEC library 
also includes post-helium core flash phases of evolution 
and has been adjusted to a revised solar abundance. In 
Figure we overlay isochrones for a 12 Gyr old popu¬ 
lation with [Fe/H] values of —1.9,—1.8, and —1.7. The 
best-fitting isochrone to the mean color of the RGB has a 
metallicity value of [Fe/H] = —1.8± 0.1; the uncertainty 
is based on the metallicities of the same age isochrones 
that span the approximate width of the RGB. The —1.8 
value seems reasonable as it is lower than the present day 
gas-phase oxygen abundance measurement of —1.52 (as¬ 
suming [O/Fe] ~ 0), and consistent with stellar metallic- 
ities spectroscopi cally measured in low-mass dwarfs (e.g., 


Kirby et al.|2010 ). We experimented with using the Dart- 
mouth stellar evolution isochrones which are known to 
provide a good fit to the RGB in low-metallicity sys¬ 
tems. The best-fitting Dartmouth isochrones had a sim¬ 
ilar metallicity value and range of —1.9 ± 0.1, provid¬ 
ing a consistency check on the values from the PARSEC 
isochrones. Using an [Fe/H] value of —1.8 ±0.1 in Equa¬ 
tion [^yields Mv = 0.52 ± 0.20, corresponding to a dis¬ 
tance modulus of 26.04 ± 0.21 mag. 

As noted above, there has been significant work done to 
calibrate the luminosity of the RR Lyrae stars. For com¬ 
parison, we consider two additional calibrations from the 
literature. First , we utilize the calibration from|Clemen-| 


tini et al. (2003) based on photometry and spectroscopy 
ot more than 100 RR Lyrae stars in the Large Magellenic 
Cloud: 


Mv = (0.866 ± 0.085) ± (0.214 ± 0.047) * [Fe/H]. (9) 


Using the same [Fe/H] value of —1.8, Equation yields 
Mv — 0.48 ± 0.12, corresponding to a distance modulus 
of 26.07 + 0.13 mag. 

Seco nd, we consider the calibration from jBono et al.| 
(20031 with uses an updated theoretical prescription of 
the pulsating modes of RR Lyrae stars: 


Mv = (0.718 ± 0.072) + (0.177 ± 0.069) * [Fe/H]. (10) 

This calibration yields an absolute magnitude of Mv = 
0.40 ± 0.15, corresponding to a distance modulus of 
26.16 ± 0.15 mag. The distances from these additional 
calibrations agree within the uncertainties to the distance 
calcul ated based on the calibration from ICarretta et al.l 
(20001 above in Equation 


3.2. HB Distance Measurement 

The luminosity of the red HB stars can also be used as 
a distance indicator. Similar to the RR Lyrae stars, the 
calibration for the HB stars has been done in the Johnson 
filters. Thus, we use the same transformations described 
in Equations and 1^ to convert the ACS magnitudes to 
V and I band magnitudes. To measure the luminosity of 
the HB, we fit a parametric LF to the observed distribu¬ 
tion of stars in the magnitude and color range of the HB 
feature, without any correction for foreground extinction. 
This is a similar maximum likelihood ap proach used to 


measure the TRGB in many systems (e.g.. 

Sandage et al. 

1979 

Mendez et al.|2002 Makarov et al.|200t)||Rizzi et al. 

2007 

accoi 

1, based on a probability estimation that takes into 
ant photometric error distribution and completeness 


using artificial star tests. We assumed the following form 
for the theoretical LF: 


P = (11) 

where A, B, and C are free parameters. We selected 
stars in the HB identified by eye in the CMD to use in 
the fit. Specifically, stars in the V mag range of 25 — 28 
and within a V—1 color range of 0.2 — 0.6 were used. The 
HB luminosity was measured to be V= 26.62 + 0.01 mag, 
from which we subtracted Av = 0.071 mag to correct for 
extinction. 
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Table 4 

Distance Measurements to Leo P 


Method 


(m-M)*^ Distance (Mpc) 


Calibration 


TRGB (LBT) 

26 

zD.iy_Q 5Q 

79 + U.i4 

-^•'^-0.40 

Rizzi et al. 

(2007 


Carretta et al. 

( 

20001 

TRGB (HST) 

zD.iy_Q 5Q 

79+0.14 

-^•'^-0.40 

Rizzi et al. 

2007 


Carretta et al. 


20001 

Horizontal Branch 

26.05 ±0.20 

1.62±0.15 



Carretta et al. 

' 

20001 

RR Lyrae 

26.04±0.21 

1.61±0.16 



Carretta et al. 

' 

20001 

RR Lyrae 

26.07±0.13 

1.64±0.10 



cl 

em 

entini et al. 

' 

20031 

RR Lyrae 

26.16±0.15 

1.70±0.12 




Bono et al. 

( 

20031 


Note. — TRGB distance was measured from LBT optical imaging of the resolved stellar populations (|McQuinn et al.||2013|l and 
confirmed via the HST dat a presented here . For t he TRGB distance measured from the HST data, we adopt ttie uncertainties trom the 



agrees 


, Lyrae 


slightly ' 


uncertainty. 


We use the HB calibration from Carretta et al. (2000| 
to determine the absolute magnitude ot the MB: 


Mv = (0.13± 0.09) *([Fe/H]-b 1.5)-t (0.54±0.04) (12) 

In the same manner described for the RR Lyrae stars, we 
assume an [Fe/H] value of —1.8 ±0.1. The HB luminosity 
based on the above calibration yields My = 0.50 ± 0.20, 
in excellent agreement with the range in calibrated mag¬ 
nitudes of the ten RR Lyrae stars. Based on this HB 
luminosity, the distance modulus is 26.05 ± 0.20. 

3.3. Adopted Distance Measurement 

Table a summarizes the distance measurements to 
Leo P based on the TRGB from the LBT and HST data 
sets and RR Lyrae and HB measurements from the HST 
data. The distance moduli determined from the different 
calibrations and technique s agree within the un certain- 
ties. The calibrations from Carretta et al. (20001 provide 
consistency across the RR Lyrae, HB, arid TRGB dis¬ 
tance techniques, allowing us to compare the three meth¬ 
ods without introducing a systematic uncertainty. These 
three techniques yield distance moduli of 26.04 ± 0.21, 
26.05 ± 0.20, and 26.19 ^q;5q respectively. The consis¬ 
tency of the RR Lyrae and HB distance measurements 
reinforces the conclusion that the TRGB identified in 


the CMD from McQuinn et al. (20131, similarly to the 
HST data presented here, is indeed a false tip. Further¬ 
more, all of the distance measurements are within the 
uncertainties estimated via Monte Carlo simulatio ns of 
synthetic populations from McQuinn et al. (2013) sup¬ 


porting this approach for determining uncertainties on 
the TRGB in this very low-mass regime. 

The HB and RR Lyrae distanc es based on the calibra¬ 
tion from Carretta et al. (2000) are very similar. We 
adopt the distance modulus ot 26.05 ±0.20 mag from the 
HB stars due to the slightly lower uncertainty, yielding a 
distance of 1.62 ± 0.15 Mpc. The refined distance places 
Leo P in the NG C 3109 association of dwarf galaxies 
(|Tully et al.|2006|) incl uding Antlia, Sextans A , and Sex- 
tans B as proposed by McQuinn et al. (2013). The clos¬ 
est of these galaxies is Sextans B at a distance of ~ 0.4 
Mpc. In Tablej^ we include updated distance dependent 
parameters using this revised distance measurement. 

4. STAR FORMATION HISTORY METHODOLOGY 

The SFH was mea sured using the numerical CMD fit¬ 
ting program MATCH ( Dolphin]|2002 ). Briefly, MATCH uses 


a prescribed IMF and stellar evolutionary isochrones to 
create a series of synthetic simple stellar populations 
(SSPs) of different ages and metallicities. The synthetic 
SSPs are modeled using the photometry and recovered 
fractions of the artificial stars as primary inputs. The 
modeled CMD that best-fits the observed CMD based 
on a Poisson likelihood statistic provides the most likely 
SFH of the galaxy. 


Th e solutions are based on a Kroupa IMF (K roupa 
2001|) and the newly availa ble PARSEC models (jBres- 


san et aL]|2012 Chen et al. 2014 Tang et al.||2014p . We 


assume a binary traction of 35% with a flat secondary 

distribution. While SFH solutions have been shown to 
be fairly insensitive to the choice of binary fraction, as¬ 
sumed fraction values b etween one and two t hirds tend to 
improve the overall fit (Monelli et al.|2010a). Extinction 
is a free parameter fit by MATCH. Both foreground and 
internal extinction can broaden the features in a CMD, 
but they are expected to be low for Leo P based on its 
high Galactic latitude and low metallicity. For eground 


extinction is estima ted to be Ag = 0.086 mag (Schlafly 
& Finkbeiner 2011) compared with the best-fitting ex¬ 


tinction value of Aij’ 475 w = 0.12 ±0.05, corresponding to 
both foreground and internal extinction which is taken 
into account in the SFH solution. The gen eral agree¬ 
ment between the foreground extinction from [Schlafly fc| 
Finkbeiner (2011) and the best-fitting extinction value 
trom MATCH provides an additional consistency check on 
the SFH solution. Because the photometric depth of the 
data does not fully constrain the metallicity evolution of 
the systems, we assumed that the chemical enrichment 
history, Z{t), is a continuous, non-decreasing function 
over the lifetime of the galaxy. Uncertainties on the SFHs 
take into account both s ystematic uncer tainties from the 
stellar evolution models ( Dolphin|20 12 ) and random un- 
certainties due to the fiiiite numbe 
polphin|[M^ . 


rer o: 


stars in a CMD 


Distance is a free parameter fit by MATCH. MATCH can be 
run without distance constraints to determine the best¬ 
fitting distance to the CMD by the models. This dis¬ 
tance can be compared with the measured distance from 
an independent approach as a consistency check on the 
solutions. The final SFH solution is derived by fixing the 
distance to be the best measured value. In the case of 
Leo P, the best-fitting distance modulus by the PARSEC 
models is 26.15±0.05 which is in good agreement with the 
measured values in Table We experimented by fixing 
the distance modulus to the different values derived from 
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WFC475W-WFC814W 


WFC4.75W-WFC814W 


Figure 7. The observed and modeled CMD for Leo P. (a) Ob¬ 
served Hess diagram of Leo P from photometry of the HST images, 
(b) Modeled Hess diagram based on the best-fitting SFH solution 
to the observations with the PARSEC library. For both (a) and (b) 
the grayscale is based on the number of stars in each bin. (c) Resid¬ 
ual Hess diagram between observed and modeled CMD of Leo P. 
(d) Residual significance between the observed and modeled CMD 
of Leo P. The grayscale spans 5(t and is based on the significance of 
each bin in the residual (panel c) relative to the standard deviation 
of a Poisson distribution. The final panel (d) highlights the good 
agreement between the modeled CMD and observed CMD. 


the HB and RR Lyrae stars. These tests showed that the 
modeled CMD with a smaller adopted distance modulus 
of 26.05 is an improvement over a larger distance mod¬ 
ulus of 26.19. The final SFH solution for each model 
is based on fixing the distance modulus to the adopted 
value of 26.05. 

In Figure 0 we present the observed CMD, modeled 
CMD using the PARSEC isochrones, residual CMD, and 
residual significance from MATCH. The modeled CMD is 
an excellent ht to the observed CMD of Leo P, with a 
corresponding value of 1.26. The residuals are most 
notable in some of the typical areas including the width 
of the RGB and the RC. While small, the residuals high¬ 
light the differences between the observations and the 
model from the isochrones. The differences in the SFH 
implied by the residuals are well within the measured 
uncertainties and do not impact our overall conclusions. 

The SFH solution depends on the ability of the stellar 
evolution isochrones to accurately model the luminosity. 



good match to the observed CMD, with slightly larger 
discrepancies in the RGB and RC regions. Specifically, 
the Padova models expected RGB stars that were bluer 


than observed, while both libraries produce more HB 
stars than observed. The corresponding values for 
the Padova and BaSTI libraries were 1.34 and 1.37 re¬ 
spectively. Note that the BaSTI models do not include 
stars with ages < 30 Myr. We present the SFH results 
for all three libraries in the following section. 


5. THE STAR FORMATION HISTORY OF LEO P 

The left panel of Figure [^presents the cumulative SFH 
of Leo P based on the best-fitting modeled CMD from 
MATCH using the PARSEC models with both random and 
systematic uncertainties. The fraction of stellar mass 
formed at each epoch is listed in Table Compari¬ 
son solutions derived using the Padova and BaSTI li¬ 
braries are shown the middle panel, with constraints on 
the chemical evolution shown in the right panel. Given 
the range in the three solutions, the SEH is consistent 
with a relatively constant level of star-format ion activity 
over cos mic timescales as seen in more massive dirrs (e.g., 
IC 1613; Skillman et al.|2014 1. The best-fitting solutions 
for Leo P also suggest that star formation may have been 
damped after early epochs (i.e., post-reionization), with 
little SF occurring between 8—12 Gyr ago. As seen in 
the middle panel of Figure]^ the solutions from the three 
stellar libraries provide a range for the degree and dura¬ 
tion of damping. Although different in amplitude, this 
is reminiscent of the delayed star form ation seen in the 
SFH of the gas- rich LG dwarfs L eo A ( Gole et M^|2007 l 
and DDO 210 (Cole et al. 2014). Regardless, from the 
SFH it is clear that Leo P 3xperienced star formation at 
t > 10 Gyr, in agreement with the identification of ten 
RR Lyrae stars. Over the last ^ 4 Gyr, the solutions 
from all three libraries show continual SF at an approx¬ 
imately constant rate. 

The stellar mass of Leo P can be estimated by inte¬ 
grating the S FH over time and ass uming a 30% recy¬ 
cling fraction (Kennicutt et al. 1994). We use the PAR¬ 
SEC solution for this calculation and find the present 
day stellar mass in Leo P to be 5.6l5;g x 10® Mq, re¬ 
ported in Table The resulting stellar mass-to-light ra- 
tio is Ma,/Lv = 1.25, based on a solar My = 4.86 ± 0.02 
(Pecaut & Mamajek 2013). The stellar mass derived 
from the SH'H is in good agreement with the previous 
estimate of 5.7('(i'8 x 10® Mq from ground-based LBT 
images ( McQuinn et al.[2013) which u sed the M/L ratio 


formalism from Bell &: de Jong 
Regardless of the model, the 


( 2001 ). 
dest-htting chemical evo¬ 


lution histories show very little enrichment in the stars 
with [M/H] values ranging from —1.7 to —1.6. These val¬ 
ues are close to our estimate of [Ee/H] of —1.8 ±0.1 from 
PARSEC isochrone fits for the oldest stars used in the 
distance calibrations of the HB and RR Lyrae stars. Rel¬ 
atively slow chemical evolution has been noted in other 
low-mass galaxies such as Leo A, IC 1613, and DDO 210 


(Cole et al.|2007 Skillman et al.|2014 Cole et al.|2014 ) 
but Leo P snows less chemical evolution than any of these 
systems. A follow-up investigation of the chemical evolu¬ 
tion history of Leo P will be presented in a future paper 
(K. B. W. McQuinn et al. in preparation). 


6. LEO P AS A PROBE OF EVOLUTIONARY SCENARIOS 
6.1. Comparison Sample and Present-Day Luminosities 
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Figure 8. Left vane l: Best-fitting c umulative SFHs for Leo P derived using the PARSEC stellar evoluti on libraries. Ra ndom uncertainties 
are plotted in yellow (|Dolphin|2013|l; combined random and systematic uncertainties are plotted in grey (|Dolphin|2012||. Middle panel: For 
comparison with the l-'AKbiiiO moaels, we present the best-fitting cumulative SFHs using two additional stellar libraries: Padova libraries 
(blue), and BaSTI libraries (red ). Random un certainties are plotted in shaded colors for each library. The systematic uncertainties derived 
for the PARSEC libraries using ||Dolphin|2012|| (shown in the left panel) are supported by this comparison. Right panel: Chemical evolution 
[M/H] solution for Leo P including random uncertainties. For clarity, we plot limited modeled values of stellar metallicity during periods 
of star formation. The chemical enrichment solutions show good agreement ranging from [M/H] of —1.6 to —1.7 across all ages. 


Table 5 

Cumulative SFH of Leo P 


Fraction of Stellar Mass 
Formed by each Epoch 


fio.io 


flO.05 


fio.oo 


fg.gs 

n 1 O+0-47 
U.l»_o_oo 

fg.go 

n 47+0-25 
-0.18 

fe.so 

-L-UU_o_oo 

fe.TS 

-*-•0^-0.00 

fe.TO 

-L-UU_o_oo 

fe.es 

-L-UU_o_oo 

fe.eo 

-L-UU o oo 


Note. — The fraction of stellar mass formed prior to each log 
time bin based on the best-fitting SFH using the PARSEC stellar 
evolution library. Uncertainties include both random and system¬ 
atic uncertainties. The SFHs are derived assuming the fraction of 
stellar mass formed is zero at log(t) = 10.15 and unity at log(t) 
= 6.6. Integrating the SFH over all time bins, the total stellar 
mass formed is 8.6 x 10^ Mq. Assuming a 30% recycling fraction, 
the present day stellar mass in Leo P is 5.6i^^ g x 10^ Mq. The full 
cumulative SFH is published in its entirety in the electronic edition 
of this journal. A portion is shown here for guidance regarding its 
form and content. 

To provide context for Leo P’s properties, we compile 
a comparison sample of very low-mass, low-luminosity 
galaxies in the nearby universe, listed in Table[^ By both 
necessity and design, the sample is heterogeneous, in¬ 
cluding both gas-poor dSphs satellites of the Milky Way 
and gas-rich dirrs that lie outside the virial radius of the 


Milky Way but still inside the LG zero velocity bound¬ 
ary. On the one hand, there are very few known gas-rich, 
star-forming galaxies in this mass range, despite predic¬ 
tions that they should be the most common galaxy struc - 
ture in the nearby universe (e.g., Haynes et al. 2011). 
As noted in the Introduction, these galaxies have proven 
elusive observationally. Thus, to build a larger sample 
of galaxies in this mass regime, we must include some 
of the gas-poor low-mass dSphs detected in close prox¬ 
imity to the Milky Way. On the other hand, the dSphs 
present very different environment-driven histories than 
Leo P and other dirrs. Thus, these galaxies provide a 
contrast to evolution in isolation allowing us to explore 
possible environmental differences on evolution in this 
mass regime. 

As an extension to our comparison, we also include 
three dSphs with even lower-lumnosities. The discov¬ 
ery and study of a population of very low-mass Milky 
Way satellite dSph galaxies, colloquially called “ultra¬ 
faint dwarf” galaxies, have extended the constraints of 
structure formation and evolution in the universe. These 


galax ies are very low-luminosity {My > —8; Martin et al. 


metal-poor, stellar population (e.g., |Sand et al. 2009 

20101 Okamoto et al. 2008 |2012t de Jong et al.||2008a 

Hughes et al.|2008 

IMartin et al.|2008l I'rebel et al.|2010 

Morris et al.||2010| 

:lirbv et al. 120101 1201 la| IBrown et al. 

2014| Weisz et al. 

2014ap. The mass-to- 

ight ratios of 
matter domi- 

these systems mdic 
nated (i.e., M/L > 

;afe that they are dark 

100 Kleyna et al. 112005 

Munoz et al. 

2006 Martin et al 

2007 Igimon &: Geha 

2M7), setting 

them apart from globular clusters of simil 

ar luminosity. 


Discovered from slight over-densities in the SDSS data 
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Table 6 

Comparison of Leo P with Nearby Very Low-Mass Galaxies 


Galaxy 

Mv 

(mag) 

Dist. 

(kpc) (10® 

M* 

Mo) 

Ref. 

Gas-Rich Low-Mass Dwarfs 


Leo A 

-12.1 

798 

60 

1 

DDO 210 

-10.6 

977 

16 

2 

Leo P 

-9.27 

1640 

5.6 


Leo T 

-8.0 

417 

1.4 

1,3 


Low-Mass dSphs 



Draco 

-8.8 

76 

2.9 

1,4 

Ursa Minor 

-8.8 

76 

2.9 

1,5 

Canes Venatici I —8.6 

218 

2.3 

1,4 


Very Low-Mass dSphs 



Hercules 

-6.6 

132 

0.37 

1,4 

Leo IV 

-5.0 

154 

0.19 

1,6 

Canes Venatici II —4.9 

160 

0.079 

1,4 


Note. — All values for Leo P are from this work. Stellar Masses 
for the comparison sample are based on a stellar mass-to-light ratio 
of unity. 


References: 

IMcConnachie (1201211: (2) ICole et al.ldi 

de Jone et al.l i 

2UU8b||: (4) liViartm et al. 

]2UU8J; [hj 

tlatzidimitriou| (|1 

995|l; (Ojlde Jong et al.|l|20 



S im 

in &| 


set, it appears that these very low surface bright ness 


galaxies are an extension of dSphs to lowe r masses (Be- 
2007 Clementini et al.||2012 ). 


lokurov et al. 


from Table ^ the present-day luminosities of the dirrs 
and dSphs span close to four magnitudes; this range in¬ 
creases to seven magnitudes when including the very low- 
mass dSphs. The dSphs are all located within ~ 200 kpc 
of the Milky Way, well inside the virial radius. The three 
comparison dir rs reside outside of t his boundary, but still 
inside the LG ( McConnachie|[20l2 ). 


6.2. Luminosity Evolution of Very Low-Mass Galaxies 

Another way to compare these low-mass galaxies is to 
consider the evolution of their luminosities with time 


based on their SFHs (e.g., Weisz et al. 2014bI. This 
can allow one to compare the initial properties of the 
sample and study whether subsequent changes are cor¬ 
related with environment. Such a comparison can be 
made by converting the SFHs into luminosity evolution 
profiles. The SFHs from resolved stellar populations of 
various low-mass gal axies have been published by a num¬ 
ber of authors (e.g., Cole et al.||200^ Kirby et al.||2011a[ 


Clementini et ai.|T2Ul2 

2014al ICole et al.boid 


jxiroy et ai.yzuiia 

|l^0l4|WeTs" etITI 


[Brown et al. ] 

among others). 

Because our analysis requires tabulated results of the 
fraction of stars as a function of time, we utilize the re¬ 
sults on eight galaxies in our com parison sampl e derive d 
from HST optical imaging from Weisz et al. (2014a), 
and similarly obtained results on one galaxy (1)1 X1 210) 
from [Cole et al.| ( |2014[ ) . The SFH solutions from jWeisz 
( |2014a| )' were derived using the Padova isochrones 


et al. 


and the solution for DDO 210 from Cole et al. (2014) 
was derived using the PARSEC isochrones. However, the 
CMDs in both of these previous studies have photomet¬ 
ric depths reaching below the old main sequence turn-off 
(oMSTO), which provides robust constraints on the an¬ 
cient star formation. SFH solutions derived from differ¬ 
ent stellar evolution libraries using data with photometry 
extending below the oMST O have shown excellent agree- 
ment with each other (e.g., Gallart et al. 2005 Monelli 


et al. 2010a|b Hidalgo et al. 2011 Weisz et al. 2011). 


Thus, the systematic uncertainties between members of 
our comparison sample should be small. For complete¬ 
ness in our qualitative comparison, we also include the 
luminosity evolution of Leo P based on the SFH from the 
Padova, PARSEC, and BaSTI models, which provides a 
sense of the systematic uncertainties. 

The SFHs were used as in put to the stellar popula- 
tion synthesis technique from Bruzual & Chariot (20031 
which builds a synthetic galaxy spectrum as a function 
of time, assuming a co nstant metallicity and an IMF. 
Because the SFHs from |Weisz et al. (2014a I and Cole 
et al. (2014) were derived from H S'l' data which do not 


encompass the full optical disks of the galaxies, a correc¬ 
tion is necessary. While the SFHs are generally assumed 
to be representative of the galaxy as a whole, the stellar 
mass only traces the stars (and luminosity) in the field 
of view. As a result, the modeled present-day luminosi¬ 
ties from the SFHs are less than the measured luminosi¬ 
ties. The exception to this is the luminosity for Leo P 
whose stellar disk is contained within the observational 
footprint. To account for the disparate fields of view be¬ 
tween SDSS imaging used to determine the luminosities 
and da ta u sed to measure th e SFHs from IWeisz et al.l 
( 2014a[) and Cole et al. (2007), we normalized the results 
from |Bruzual & Chariot 1 2003) by the present-day mea- 
sured luminosity. I'his scaled the modeled luminosities 
while preserving relative changes as a function of time. 

In the left panel of Figure we present a comparison 
of the luminosity evolution (Mv{t)) for Leo P based on 
the three stellar evolution libraries. Because the best¬ 
fitting solution from the Padova models has 50% of the 
stellar mass formed in the earliest epochs, compared to 
10% — 20% from the BaSTI and PARSEC libraries re¬ 
spectively, the Padova luminosity profile also shows the 
highest luminosity at earliest times. Regardless of the 
absolute magnitude, all three solutions show a peak in 
brightness at early times followed by a slow fading in 
luminosity. 

The right panel of Figure presents a comparison of 
Mv{t) for Leo P using the PARSEC models with six 
gas-poor dSphs and three gas-rich dirrs. We chose to 
highlight the Leo P solution using the PARSEC mod¬ 
els which provided the best-fitting modeled CMD to the 
data, but our overall interpretation is valid for the pro¬ 
files derived from all three stellar libraries. Two of the 
very low-mass dSphs (CVn H and Leo IV) lie 2 — 4 mag¬ 
nitudes below the rest of the sample in both initial and 
present-day luminosity. Excluding these two systems, 
the rest of the sample shows a much narrower range in 
their initial luminosities than their present day luminosi¬ 
ties (i.e., AMy ~ 2 versus AMy ~ 5). The larger range 
at the present-epoch is mainly driven by changes in three 
galaxies. On the bright end are DDO 210 and Leo A 
which show dramatic increa s es in star formation in the 
last 6 Gyr (Cole et al. 2007 2014). On the faint end is 
Hercules which has formed far fewer stars at recent times 
than the initially comparable low-luminosity dirr, Leo T. 

Looking at the early evolution of the profiles in Fig¬ 
ure Leo P’s initial luminosity is similar to one very 
low-mass dSph (Hercules) and lower than three other 
dSphs (CVn I, Draco, Ursa Minor).This similarity sug¬ 
gests that Leo P may be what a dSph would look like if it 
evolved in isolation and retained its gas. One of the ar- 
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Age (Gyr) 


Age (Gyr) 


Figure 9. The luminosity evolution of Leo P and nine low-luminosity Milky Way satellites based on their SFHs and stellar population 
synthesis modeling. The luminosity of Leo P at early times is consistent with a number of the now gas-poor dSphs which have since faded 
to lower luminosities due to their lack of star formation. Leo P is analogous to a dSph evolving in isolation, highlighting the dramatic effect 
environment has on very low-mass galaxies. Note the smoother profile of DDO 210 is due to the larger time binning of the SFH from|Cole| 
|et al.| | [2014[ |. 


guments against a dirr evolving into a dSph is based on 
conservation of angular momentum; dirrs are typically 
rotat ion supported whe reas dSph are pressure supported 
(e.g.,|Grebel et al.|2003|). However, analysis of the gas ve¬ 
locity neldurTCedPsEows that the gas velocity dispersion 
approaches the rotation speed. As discussed in detail in 
Bernstein-Gooper et al. (2014), there is significant ran- 
dom motion of the gas in Leo P superposed on the bulk 
rotation signature. This agrees with previous findings of 
little evidence of rotational support in other very low- 
mass dirrs such a s Leo A and DDO 210 (Lo et alj|1993 


authors have attributed the removal of gas and metals in 
Draco, GVn I, and Ursa M inor to stellar feedback pro¬ 


cesses (Kirby et al. 
be involved to exp 


2011bI, but additional factors must 
ain the distinct retention of gas in 


Young fc Lo|1997 ). In contrast to higher-mass dirrs with 


clearly dehned rotation curves, the transformation of a 
very low-mass galaxy like Leo P to a dSph would not 
require a drastic loss of angular momentum. 

6.3. Similarities within Morphological Types at Recent 

Times 

From Figure the patterns in luminosity profiles 
within the two morphological groups are quite similar 
over the last ~ 7 Gyr. The dirrs show fairly steady or 
increasing luminosities whereas the dSphs show declining 
luminosities. This more recent evolution in luminosity is 
presumably dictated by environment, described by the 
morphology-density relation (i.e., the isolated low-mass 
galaxies which retain their gas continue forming stars 
while those in closer proximity to a massive galaxy do 
not retain their gas and begin to fade). This is indi¬ 
rect evidence that the removal of gas in the dSphs is 
driven by environment — not internal stellar feedback — 
although feedback may be necessary to oust the gas from 
the inner, denser parts of the galaxies into an enveloping 
diffus e gaseous halo (see, e.g., Mayer et al 2006 Simpson 


comparably low-mass systems such as Leo P and Leo T. 

6.4. Similarities across Morphological Types at Early 
Times and Implications 

Despite the heterogeneity and general divergence of the 
profiles between the morphological types at more recent 
times, there is a similarity in the first 6 Gyr regardless 
of present-day morphological classification. As seen in 
Figure there is a peak in luminosity at early epochs 
followedDy a decrease. As suggested by previ ous studies 


(Cole et al. 2007 2014 Brown et al. 2014), this pat¬ 


tern IS consistent with a quelling of star torniation post¬ 
reionization. In this scenario, ionizing photons f rom the 
first epochs of star forma tion {z = 11.3 ± 1.1; [Planck 
Collaboration et al. 2014), may heat the outer disks of 
the galaxies in the mass and distance ranges probed, 
slowing or halting further gas accretion and causing an 
eventual (althou gh sometimes tempo rary) decline in star 


eventual (aitnou gn sometimes tempo rary) decline m stai 
formation (e.g., [Ohorbe et al. 2015|). As noted above 


et al. 2013[ [Milosavljevic fc Bromm|201^ and references 

therein) for more swift and efficient removal. Previous 


the galaxies at closer distances to the hostile environ¬ 
ment of a massive galaxy are more likely to have the 
hot coronal gas removed by stripping mechanisms. The 
more distant, isolated galaxies retain the gas heated by 
reionization, which is likely at extended radii in the dark 
matter halo. Through atomic line cooling processes over 
Gyr-long timescales, this hot coronal gas may eventu¬ 
ally condense onto the disk and reignite star formation, 
consistent with the increases noted for the more isolated 
galaxies (Leo P, Leo A, DDO 210, Leo T). 

The impact of reionization on an individual galaxy is 
likely dependent on the distance to the local source of 
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the ionizing photons, ho wever it is unclear what distance 
scales are involved (e.g., Weinmann et al.||2007 and ref¬ 
erences therein). Excluding Leo P, the distances probed 
by our comparison sample range from 76 kpc — 1000 kpc; 
within this range, the early evolution of the luminosity 
profiles do not show differences that can be obviously 
tied to present-epoch distances. For example, contrary 
to the idea that galaxies closer to the ionizing source 
should be more impacted by reionization, the galaxies 
closest to the Milky Way (Draco, Ursa Minor, CVn II, 
Leo IV) show smaller fractional increases in luminosity 
between ~ 8 — 12 Gyr ago; the more distant galaxies do 
not, suggesting additional variables (e.g., mass) must be 
considered. The literature SFHs used to model the lu¬ 
minosity evolution have similarly small uncertainties at 
old look back times. Thus, the luminosity profiles should 
be robust. The exception, of course, is Leo P, which at 
the larger distance of 1.6 Mpc, limits the accuracy of the 
ancient SFH with the current data set. The best-fitting 
SFH suggests star formation was damped even at the 
larger distance of Leo P. However, given the uncertain¬ 
ties at older times, the SFH of Leo P does not preclude a 
fairly constant SFR as an alternative evolution for Leo P. 
If true, this would imply the damping of star formation 
occurs inside of ~ 1.5 Mpc scales for a LG-type environ¬ 
ment. 

The impact of reionization on early star formation is 
also predicted to be dependent on the mass of a galaxy. 
Theoretical simulations fine-tune a “filtering mass” be¬ 
low which reionization universally quenches star forma¬ 
tion. While the “filtering mass” g enerally falls within 
the mass range pr obed here (e.g., Gnedin 2000[ Susa 


& Umemura 2004), we do not see evidence ol univer¬ 
sal quenching, m agreement with previous results (e.g 


Grebel fc Gallaghen|2004[ IMonelli et al.|[^10a| |Hidalgo| 
et al. 2011 Weisz et al. 2015p . Yet, observational stud¬ 


ies have reported something akin to this “knife-edge” 
in star formation within the temporal limitations of the 
GMD analysis in a few dSphs of even lower-masses (i.e., 


M* ^ 10^ Mq Brown et al. 


2014). There is still con¬ 


tinuing debate on whether the cessation of star forma¬ 
tion in the aforementioned systems can be sol ely at¬ 


tributed to photoevaporization from reionization (Brown 


et al. 2014) as these galaxies are also susceptible to 


environment ally-driven removal o f gas, and possibly stel¬ 
lar feedback ( Kirby et al.||2011b |. 

Finally, we note that while Leo P has initial luminosity 
properties similar to the dSphs Hercules, Draco, GVn I, 
and Ursa Minor, the overall evolution of Leo P, based on 
the best-fitting SFH, is most similar to Leo T (shifted to 
brighter magnitudes). Leo T lies just outside the virial 
radius of the Milky Way, but is thought to have been 


only recently accreted {tinfail < 1 Gyr; Rocha et al. 
2012P . The similarities in the luminosity prohles support 
this hypothesis as tidal stripping of gaseous material in a 
galaxy as low-mass as Leo T becomes important at ~ 1.5 
times the virial radius of the Milky Way, or ~ 450 kpc 
(Milosavljevic & Bromm||2014). Thus, it is unlikely that 
Leo 'L would have retained its gas if it had been in resi¬ 
dence around the Milky Way for a cosmically-significant 
time. Previous authors have suggested that Leo T, as 
well as Leo A, fit in the general scenario outlined above 
for Leo P where star formation is quelled by reionization 


and reignited by late phase cold g 

as accretion dCole et al. 

2007 

Ricotti 2009[ Weisz et ai.| 

2012 Glementini et al. 


2012 ). 


7. CONCLUSION 

Using HST optical imaging, we have studied the stellar 
populations and evolutionary history of Leo P. We refined 
the distance measurement to Leo P to 1.62 ± 0.15 Mpc 
based on HB stars and ten RR Lyrae candidate stars, 
identified from their light curves. Using a GMD-fitting 
technique, we have reconstructed the SFH and chemical 
evolution of Leo P. The best-fitting SFH solution favors 
star formation at the earliest epochs, followed by a period 
of quiescence, and a relatively constant star formation 
rate at recent times. 

The closest galaxies provide the most detailed con¬ 
straints on the evolution of galaxies at the faint-end of the 
LF. Using luminosity evolution profiles modeled from de¬ 
tailed SFHs, we note a similar pattern between low-mass 
dirrs and comparably low-mass dSphs at early times. 
The luminosity profiles of both morphological types sug¬ 
gest a general quelling of star formation post-reionization 
with no obvious distance dependency, after which marked 
dissimilarities arise. 

There are a number of possible explanations for the 
quelling of star formation and subsequent evolutionary 
differences in this low-mass regime, but the scenario 
which fits the observed trend includes the following gen¬ 
eral framework: (i) heating of the outer gaseous disks of 
low-mass galaxies by reionization at the earliest epochs, 

(ii) the cessation of gas accretion onto the galaxies, which 
provides a governor on the total gas mass of each system, 

(iii) in the case of dSphs, the stripping of the hot coronal 
gas in the hostile halo environment of the host galaxy 
(e.g., the Milky Way), (iv) in the case of isolated dirrs, 
the slow cooling of the ionized gas which eventually con¬ 
denses and becomes available for further star formation. 
The similar initial luminosities of low-mass galaxies stud¬ 
ied provides supportive evidence that the removal of gas 
in post-reionization epochs has a strong environmental 
dependence for galaxies in this mass regime. Addition¬ 
ally, because both isolated and non-isolated galaxies will 
suffer the loss of material from their shallow potential 
wells due to stellar feedback processes during epochs of 
star formation, stellar feedback alone cannot explain the 
present-day luminosity differences between the galaxies. 

Whether or not the SF in Leo P at early epochs was 
temporarily damped by re-ionizing photons remains an 
open question. The SFH suggests there may have been 
a cessation in activity from 9—13 Gyr ago. On the 
other hand, within the SFH uncertainties the solution is 
also consistent with a relatively constant SFR scenario 
with no quenching. If a more constant SFR is favored, 
this would imply that the impact of reionization is sig¬ 
nificantly smaller outside of a ~ 1.5 Mpc radius from a 
LG-type environment. Despite this ambiguity, it is clear 
that Leo P contains a population of old stars and has also 
retained a significant gas reservoir and continues making 
stars at a fairly constant rate. The initial luminosity of 
Leo P is similar to dSphs, suggesting that Leo P is what 
a low-luminosity dSph would look like if evolved in iso¬ 
lation without losing its gas and its ability to form stars. 
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